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Abstract 

Currently, certification of engine designs for resistance to bird strike is reliant on physical tests. Predictive 
modeling of engine structural damage has mostly been limited to evaluation of individual forward section 
components, such as fan blades within a fixed frame of reference, to direct impact with a bird. Such 
models must be extended to include interactions among engine components under operating conditions to 
evaluate the full extent of engine damage. This paper presents the results of a study aim to develop a 
methodology for evaluating bird strike damage in advanced propulsion systems incorporating hybrid 
composite/metal structures. The initial degradation and failure of individual fan blades struck by a bird 
were investigated. Subsequent damage to other fan blades and engine components due to resultant violent 
fan assembly vibrations and fragmentation was further evaluated. Various modeling parameters for the 
bird and engine components were investigated to determine guidelines for accurately capturing initial 
damage and progressive failure of engine components. Then, a novel hybrid structure modeling approach 
was investigated and incorporated into the crashworthiness methodology. Such a tool is invaluable to the 
process of design, development, and certification of future advanced propulsion systems. 

1 Introduction 

The threat of foreign object impact (FOI), such as bird strike or hail strike, on turbofan engines is a major 
concern for commercial airlines and engine manufacturers. The FAA has set specific airworthiness 
standards under Federal Aviation Regulation (FAR) 33.76 requiring that engines be capable of 
withstanding impact with birds ranging from 0.8 to 8 lbs. without sustaining damage that poses a fatal 
threat to passengers and crew. To meet these standards, jet engine manufacturers rely heavily on a 
multitude of costly physical tests requiring the destruction of full-scale engines by simulating bird strikes 
using birds or bird surrogates in accordance with FAA guidelines. 

Bird strike on a turbofan blade can be defined as a complex, transient fluid-solid interaction between the 
soft bird material and the fan blade [1]. Upon impact, the internal solid structures of the bird are rapidly 
deteriorated as the resultant shockwave propagates through its body, breaking internal bonds of the bird 
and transforming its material properties to that of a viscous fluid contained within a shelled control 
volume. This causes pressure to be built up within the bird material and, as the shockwave progresses, 
leads to the formation of pressure release waves that force the transient fluidic matter out of the initial 
control volume enclosure [2]. The bird material expands nonlinearly beyond the initial point of impact, 
creating a larger contact region and causing greater potential for damage to additional fan blades as well 
as other engine components outside of the initial impact zone. 

Several types of damage occur in the forward fan section of an engine once subjected to an FOI. One type 
of damage is that resulting from a bird directly impacting several fan blades. Disintegration and 
degradation of the fan blades can cause decreased engine performance and further damage to other 
components as fragments or detached blades are ingested into the engine. 

Another significant source of damage is that resulting from contact between the fan blades and engine 
casing leading to high tangential friction and normal impact forces along the blade tips and the inside wall 
of the engine casing. Minute disturbances in the fan assembly position will bring the blades into contact 
with the casing, producing high blade tip stresses as they shear against sections of the engine casing. Bird 
impact along the length of a fan blade induces a moment about the engine axel and can remove mass from 
portions of the fan assembly, shifting the center of mass away from the axis of rotation. This can induce 
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vibrations that move the fan blades through the small blade tip toleranee. Fragments resulting from blade 
and easing damage are ingested into downstream seetions of the engine and ean eause further destruetion. 

Capturing the deterioration of engine eomponents resulting from violent eontaet between the fan blade 
tips and the engine easing is essential for eharaeterizing damage in turbofan engines. This is a partieularly 
ehallenging task for modern fan blades utilizing eomposites with reinforeing metal eomponents, with the 
resulting part ealled a hybrid strueture. Tools for aeeurately predieting damage in these struetures are not 
well established. New methods are needed to faeilitate the design and analysis of future propulsion 
designs ineorporating these advaneed struetures. 

The goal of the eurrent study is to develop a methodology for predieting bird strike damage in jet engine 
forward fan seetions ineorporating the true nature of struetural dynamies and destruetive interaetions. 
Critieal features under investigation are the degradation and disintegration of the turbofan blades and the 
engine easing. Three-dimensional explieit finite element analysis models will be employed to simulate the 
highly nonlinear, transient response of the forward engine seetion resulting from impaet and eapture the 
assoeiated turbofan blade damage, fan assembly vibrations, and subsequent engine easing damage. A 
methodology for predieting these damage effeets is imperative to faeilitating the design and eertifieation 
proeesses of future advaneed propulsion systems. 

2 Problem Approach 

A review of the literature shows several approaehes employed through previous researeh programs to 
prediet bird strike damage to jet engine fan assemblies. Several researehers have eondueted in-depth 
studies of damage to an individual fan blade resulting from bird impaet to assess the response and damage 
of the strueture [3-5]. In addition, some investigators have studied the effeets of impaet between a bird 
and multiple (three or four) blade eonfigurations in aeeordanee with the distributed impaet that oeeurs 
when a bird strikes several blades as it passes through the front seetion of the engine [6, 7]. Lastly, there 
are few models that inelude the entire fan assembly arranged about a eentral hub in whieh the fan blades 
rotate about a rigid point. 

While these models eapture the damage that oeeurs to a blade or blades direetly impaeted by a bird, they 
do not ineorporate the engine easing or a deformable axel into the simulation and the highly interaetive 
role that their dynamie responses play in sustaining or eontaining damage. Therefore, these models are 
unable to eapture the post-impaet damage responses that may result from abrasion of the fan blade tips 
onto the engine easing. Also, these studies are often limited to individual eases in whieh one simulation is 
developed to mateh experimental or analytieal data without eonsideration for the effeet that modeling 
ehanges may impose on predieted results. 

The approaeh ehosen for this study was to model the bird impaet damage for a jet engine forward seetion 
eapturing several key features: 

Initial impaet between the bird and fan blades and the resulting total damage. 

Mass unbalaneing, leading to fan assembly and main axel vibration and post-impaet dynamie 
damage. 

Blade tip shearing against the engine easing and the resulting damage both the blade and the 
easing. 
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Additionally, it is aimed to explore modeling approaehes that affeet the predieted initial and final damage 
in order to establish general guidelines for bird strike damage assessment on engine forward seetions and 
any possible eonelusions eontributing to the design or eertifieation proeedures. 

3 Problem Complexities 

The diffieulties of modeling bird strike damage on a full fan seetion are numerous even for a simplified 
model. In addition to the eommon ehallenge of aeeurately modeling progressive material failure, the 
following major eomplexities must also be addressed: 

3.1 Bird/Fluid Modeling 

An important eomplexity was modeling the bird aeeurately using an approaeh that ean eapture large 
deformations and separation of the fiuidie bird material following initial impaet without ineurring heavy 
element distortions that ean lead to mesh entanglement and ultimately divergenee of the analysis. To 
avoid this, a meshless formulation was utilized for the bird. This approaeh is diseussed in detail in seetion 
5. 


3.2 Contact Modeling 

Essential interaetions required to be ineorporated into models ineluded bird-blade, bird-easing, blade- 
easing, blade -blade, and blade-self eontaet. Contaet between the bird and engine eomponents had to 
aeeurately eapture the interaetion forees between the bird and eomponent surfaees over time. Sinee 
eontaet between interaeting parts due to material failure in element eliminated areas was lost, eontaet 
regions had to be redefined over time as the simulation progressed. FE models had to be able to 
adequately impose eontaet between eroding parts as eontaet interfaees eontinued to evolve. 

In addition, further eomplexity arose due to the elose proximity and high veloeity of the blade tips relative 
to the easing elements. In post impaet, blade elements penetrated the easing due to oseillations assoeiated 
with unbalaneing in the fan assembly. Contaet models were utilized that deteeted this penetration and 
applied subsequent interaetion forees to remove the penetration and thus model the eontaet between 
elements. However, if the time step was too large, the blade elements penetrated the easing beyond the 
limits that the eontaet formulation was eapable of deteeting and failed to impose the eonstraint, eausing 
unrealistie penetrations between the parts. To prevent this, the initial time step was lowered until it was 
observed that exeessive element penetration no longer oeeurred for the range of fan assembly veloeities 
under investigation. 

3.3 Fan Unbalancing and Vibrations 

Impact between the bird and fan blades at a distance away from the center of the hub moved the fan 
assembly off center, causing it to oscillate about the engine’s axis of rotation. In addition, if the blades 
fractured, the center of mass shifted permanently, triggering a strong rattling effect as the fan blade tips 
hit the casing and caused damage to both the blades and the casing. Suitable models were required to 
capture damage not only from direct impact, but resulting from post-impact damage to the rest of the 
structure due to unbalancing and vibrations. 
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3.4 Pre-stressing 


The angular velocity of a jet engine fan section causes centripetal forces along the length of fan blades 
that pull them toward the center of the hub and maintain the circular path of the fan stage around the axis 
of rotation. Stresses are greatest at the blade root and decrease along the length of the blade until reaching 
zero at the tip. In order to accurately predict damage, it was necessary to assess the state of stress and 
strain in the blades as a result of both the bird impact and the initial state of stress due to centripetal 
acceleration. 


4 Soft Impact Theory 

At velocities around and above 50 m/s, a bird impacting an object with a much higher stiffness induces a 
shockwave that propagates in the direction of relative velocity through its body and breaks its internal 
bonds. This leads to a multi-stage problem as the bird’s properties, past the shock, change from solid to 
fluid. The several facets of the impact event can be characterized by the following main stages: (1) solid- 
structure impact, (2) solid-fluid transition, (3) fluid-structure interaction, (4) nonlinear growth of impact 
region, and (5) structural damage. 


Originally, the bird mass is solid, and the initial impact with a fan blade can be labeled as a “hard” impact 
in which contact forces are highly localized about a relatively small impact region. The velocity of 
material over the contact area instantly decreases, and a shock pressure develops in both the projectile and 
the target as material is compressed. Assuming both objects behave elastically, an expression for the 
shock pressure, or Hugoniot pressure, P_H, induced in the bird can be derived from: 






yp’u:+p^u. 


( 1 ) 


where Pq and p* are the densities of the bird and the blade, respectively, and are the shockwave 
speed in the bird and blade, and J/q is the relative impact velocity [1]. In the case of bird strike, the 
following will be valid: p*U^ » P(jf^j, and subsequently the above equation can be simplified to: 








( 2 ) 


The shock pressure created inside of the bird is greater than its material strength. Therefore, the structural 
bonds of the bird are broken, and the solid material undergoes a transition to a fluid. 


Concentration of the shock pressure along an interface, bound by the speed of sound in the object, 
produces a shockwave. The shockwave then propagates at a speed much faster than the relative bird/blade 
velocity, at which the bird material properties are rapidly changed into fluid. Consequently, a large 
pressure gradient is formed across the boundary of the surrounding air and the control volume containing 
the bird. Pressure release waves force the fluidic matter away from the bird center, while the shock 
pressure decays to a lower steady pressure [1]. This steady pressure, is defined by: 

( 3 ) 
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The theoretical duration of the impact event is the time it takes for the end of the bird to reach the target. 
Since the post shock gelatin like material can be considered incompressible, its speed does not decrease 
significantly throughout the remainder of the impact. Therefore, the initial velocity may be used to 
estimate the impact window as: 


T = L/U^ (4) 

where L is the length of the bird and is the initial bird velocity. 

In high velocity bird or porous Newtonian fluid-compliant gelatin impacts (Fig. 1), a sharp rise in the 
internal pressure of the soft body can be detected as the projectile first comes into contact with the target. 
This triggers a progressive change in the impactor material from solid to fluid followed by shock pressure 
release waves occurring at the solid-fluid interface. 
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Fig. 1. Pressure in a porous gelatin projectile high velocity impact [8]. 


Subsequently, a sudden pressure drop is experienced in the forward and transient fluidic portion of the 
projectile. The residual pressure then continues to gradually decrease through the remainder of the impact 
event until it finally reaches zero [1]. 

The force exerted on a target over time by a soft impactor such as a bird is directly related to this time- 
varying pressure. The applied force can be derived using the relation pressure = force/area. Therefore, 
the peak force and steady pressure are: 


= PrA (5) 

Ps = ( 6 ) 

The force therefore behaves similar to the pressure by first rising quickly to a maximum value in 
correspondence with the initial impact, decreasing rapidly as the initial shock in the projectile releases, 
becoming steady, and then declining towards zero at the end of the impact event. A linear approximation 
of this process is shown in Fig. 2. 
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Fig. 2. Linear approximation of the force time history of a soft impact. 

The excitation force, F(t) , can be approximated using a piecewise linear function that defines the active 
force over each unique time domain. Using the Heaviside function, H(t), F(t) can be defined as 

■m 

F(t) = ^ [H(t - - H(t - tJ]F^(t) (7) 

n = l 

where m is the last time piece. The equation F(t) determines the applied force with time as different 
force equations are turned “on” and “off’ by the time -dependent Heaviside functions. Therefore, for 
t > the applied force is zero after the impact has completed. This simplified model can be 

expanded to include additional time domains describing the applied force in more detail than what is 
shown in Fig. 2 in order to more accurately capture the force-time history on the target. An application of 
this concept to bird impact with equations detailing the force over each time domain is discussed 

in [9]. 

Using this excitation force, the displacement response x{t) of an impacted target can be defined. For 
example, if the target was represented as an equivalent one -dimensional mass-spring-damper system, the 
displacement of a mass subjected to the force F(t) appeared in Eq. (7) can be defined using the 
convolution integral approach: 


x(t) = 



1 




mco^ 


sin (iijjCt — z)\dT 


( 8 ) 


where m is the mass of the target, and are the damped and un-damped natural frequencies of the 
system, respectively, and C is the damping ratio. The variable t in the equation is the integration variable. 

5 Meshless Lagrangian Bird Modeling Approach 

Several meshless finite element modeling approaches have been developed that represent objects as a 
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finite number of diserete partieles that interaet via an approximation or kernel fiinetion rather than a 
struetured mesh. The range of influenee and interaetion with other partieles ean be defined through using 
this funetion. In meshless Lagrangian approaehes, eaeh node represents a partiele with physieal properties 
sueh as mass, volume, and pressure. Eaeh partiele experienees fluid forees based on the mass and number 
of surrounding partieles and the kernel defining their interaetion [10]. Closer partieles exert more 
influenee than those that are farther away. In this manner, partiele interaetion ean be visualized as an 
inter-moleeular attraetion type foree holding the partieles together. The behavior of the partieles is 
governed by the Navier-Stokes equations, numerieally diseretized for use with the partiele approaeh [9]. 
This might negleet the possible non-Newtonian nature of the transient fluid produeed as a result of the 
fluid-solid interaetion. 

5.1 Bird Model Details 

The bird was modeled as diserete partieles using a Meshless Lagrangian approaeh on the explieit finite 
element platform Ls-Dyna. Several representative bird geometries have been proposed for both impaet 
experiments utilizing bird surrogates and finite element models of high veloeity bird impaet. One of the 
most eommon eonfigurations is a eylinder with rounded hemispherieal ends [11]. 

A similar geometry was adapted for the employed partiele modeling approaeh in thus work; the partieles 
were arranged, evenly spaeed, throughout a rounded-eylindrieal volume to approximate the shape of a 
bird body. The ratio of total length to diameter for the bird geometry was 2. 

It was neeessary to define both a suitable material model and a pressure -volume state equation for the 
bird. As diseussed above, the bird material begins ehanging from a solid to a fluid at the point of impaet 
with a blade. Material at the blade/bird interfaee is rapidly liquefied, and the resulting impaet shoekwave 
propagates through the bird mueh faster than the relative bird/blade veloeity. For the majority of the 
interaetion, the blade interaets with the gelatin-like fluid, and only interaets with the solid bird for a brief 
instant at the moment of impaet. Therefore, it is justified to model the bird using a fluid material model 
for the entire impaet event. 

An appropriate fluid material model was identified that allowed the mass density and the fluid pressure 
eutoff to be defined. In addition, a polynomial equation of state (EOS) was used that defined the inverse 
relationship between the partiele volume and pressure, given as: 

F = 



( 9 ) 

C^ = K 
= (2fe- i)q 
Cg = (fe- I)(3fe-I)q. 

The eonstant is the initial density defined in the material model, p is the instantaneous density, and K 
is the bulk modulus of the bird material. Through the definition of density, the pressure -volume 
relationship within the above EOS formulation ean be established. The bulk modulus of water (2.2 GPa) 
was used in the EOS. The initial material density used was 945 kg/m^ with a total bird mass of 1.8 kg, or 
4 lb. These properties, along with the stated length-to-diameter ratio, defined the neeessary bird geometry. 
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The variable fe is a non-dimensional material eonstant ealled the Hugoniot number or eompressibility 
eonstant. This is adopted from the elassieal linear shoekwave speed equation: 

^^0 ( 10 ) 

where is the speed of sound in the material, k is the Hugoniot number, and and are the 
shoekwave speed and initial veloeity respeetively, as defined previously. Prior studies have found k to 
equal approximately 2 for biologieal materials ineluding birds [12]. 

Simulations of the bird impaeting a rigid plate were eondueted and the pressure at the eenter of the impaet 
region was reeorded throughout the duration of the impaet. The number of partieles used to represent the 
bird was varied and the pressure -time history at the eenter of impaet was eompared. This was done in 
order to verify that a suffieient partiele density was used sueh that the bird behavior was independent of 
the number of partieles used in the simulation. 

5.2 Impact Results 

The pressure -time history of the bird impaeting a rigid plate is shown in Fig. 3. The reeorded pressure was 





Nori-dimerisionalized Time, 

Fig. 3. Pressure-time history of the bird/rigid plate impact scenario for increasingly more particles used to 

represent the bird. 
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non-dimensionalized by dividing the results by the theoretieal steady pressure given by Eq. (3). Time was 
also non-dimensionalized by dividing eaeh time inerement by the theoretieal impaet duration defined by 
Eq. (4). 

Eaeh graph was eharaeterized by a short duration pressure rise, followed by an approximate steady 
pressure state, and a final region where the pressure deeays to zero. For Fig. 6 (a) and (b), the predieted 
pressure response ineluded two peaks before settling to a steady pressure. As the number of partieles was 
inereased, the response ehanged sueh that the initial high pressure region was eharaeterized by a single 
pressure rise rather than the two initially observed. 

At around 30,000 to 35,000 partieles, the pressure -time history began reaehing a eonsistent response. 
Further inereases in the number of partieles lead to minor refinements of the response with redueed 
oseillations with the overall response ehanging only slightly. The peak pressure value was also quite 
eonsistent with an approximate non-dimensionalized pressure of 8.5. This response was repeated in 
further analyses, up to 150,000 partieles. It was eoneluded that for numbers of partieles around 35,000, 
the bird model exhibited an aeeeptable level of independenee from the partiele density. To provide 
aeeuraey and minimize eomputational time, a bird model eomprising of approximately 35,000 partieles 
was used for impaet damage analysis onto the engine forward seetion. 

6 Forward Section Impact Modeling 

The fan blades were modeled with shell elements using a general eommereial aviation fan blade 
geometry. The width at the base was 0.2 m, representing a variety of wide-ehord fan blade designs. A 
total of 36 blades were modeled with their bases fixed to a eentral hub via rigid eonneetors eonstraining 
both translational and rotational motions against the hub. A total shell thiekness of 3.3 x 10'^ m was 
assigned to eaeh blade element. An isotropie plastieity material model was used to define titanium 
material properties with yield and failure strain ineorporated into the material behavior. The pre -yield 
elastie and post-yield tangent moduli of 112.5 and 1.6 GPa were respeetively seleeted. The maximum 
strain at whieh failure would oeeur was set to 0.1, determined from typieal stress-strain eurve data for an 
aerospaee grade titanium alloy. 

Post-impaet unbalaneing and vibration were determined through assessing the deformation imposed on 
the engine axel. Translational motion along its axis was restrieted via boundary eonditions on the end 
opposite of the fan assembly. Supports were plaeed along the axel to prevent rigid body movement in the 
plane perpendieular to its longitudinal axis. These supports were modeled as rigid “sleeves” that were free 
to rotate about any axis, pivoting on a eentrally loeated referenee node, however restrieted from rigid 
body translation. This suffieiently held the axel in plaee while allowing it to bend in aeeordanee with 
forees indueed by mass unbalaneing at the fan seetion. 

The easing was modeled as an open eylinder using shell elements surrounding the fan seetion. Care was 
taken to fit the diameter of the easing elose to the fan blade tips in aeeordanee with the tip toleranees that 
are used in modern engine designs. It was essential to allow a small toleranee for the fan blades to deform 
radially as a result of pre-stressing in order to avoid initial penetrations between the blades and the easing. 
An isotropie plastieity material model was used to inelude strain rate effeets, yielding, and failure. 
Aerospaee grade aluminum alloy was used for the property definition. The finite element model as a 
whole with the easing is shown in Fig. 4. 
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Fig. 4. Explicit finite element impact environment. 


An explicit dynamic relaxation method was used to pre-stress the fan blades. This approach allowed 
stresses and displacements to be initialized using prescribed centrifugal body loads by damping nodal 
velocities by a prescribed scaling factor at each time step. The peak distortional kinetic energy was 
calculated at the first time step, and the current distortional kinetic energy was determined at each 
subsequent time step [13]. When the ratio of the current distortional kinetic energy to that of the peak 
energy fell below a set threshold, the system was determined to be sufficiently steady, and the analysis 
was considered complete. As the threshold was lowered, the solution began to improve. 

For systems with loads due to rotation, it was found to be extremely important to set this value as low as 
possible. Finding the lowest convergence ratio required a trial and error procedure. After wide-ranging 
experimentation, a value of 1.5 x 10'^ was selected. The results of the pre-stressing are shown in Fig. 5. 



Stress (v-m) 
1.141e+09_^ 
1.027e+09 
9.131e+08 
7.990e+08 
6.849e-H)8 
5.707e+08 
4.566e-H)8 
3.424e-K)8 
2.283e+08 
1.141e+08 
O.OOOe+00 




Fig. 5. Fan assembly stresses following the pre-stressing phase. 


Subsequent to the stress and displacement initialization phase, the bird impact was introduced into the 
simulations: the initial velocity of the bird was 350 m/s, while the initial rotational velocity of the fan 
assembly was selected to be 525 rad/s, corresponding to a rotation of 5000 rpm. The bird was positioned 
close to the fan blade, just outside the penalty length, to reduce needless computation. 
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6.1 Bird Impact Results 


Initially, an impact scenario without the casing was modeled to establish the response of the fan assembly 
to a bird impact without interference from the fan containment. The bird directly impacted several blades 
and induced high stresses and strains that lead to material failure at the base of two of the blades. From 
this analysis, deformation of the axel was used to establish the mode shape of the vibration response of 
the fan assembly due to mass unbalancing. The presence of the three non-translating sleeves holding the 
axel meant that there were three fixed points representing the engine axel anchor stations within the 
propulsion system. The distributed mass between the sleeves and the large fan assembly constituted a 
system with three degrees of freedom. The mode shape diagram is given in Fig. 6. This mode shape 
corresponds with the first vibration mode of a three degree of freedom system. 



Fig. 6. Mode shape of the axel vibration. 


Consequently, a bird strike event against the full fan assembly with the casing included was simulated to 
observe the damage behavior of the system. Again, the bird first made contact with a single blade and 
began deforming heavily as a result of the interaction. 


The resulting soft body impactor was sliced with portions of it coming into direct contact with three to 
four blades. The force-time history of the impact between the bird and the fan blades is shown in Fig. 7. 
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Fig. 7. Impact force between bird and fan blades. 


The spread of the bird contributed to a larger contact area between the bird material and the back of each 
blade as the blades sliced through and then collided with clusters of particles. This larger contact area 
appeared to contribute to high deformation across the entire blade, inducing plastic strains through the 
blade and at its base. Eventually, strains at the base exceeded the designated material limit, causing 
failure. This was effectively captured by the analysis in the form of element deletion along the base of 
one or two of the impacted blades. 
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Fig. 8. Impact between the bird and fan assembly. 


Following the initial impact, several events occurred that contributed to extensive damage to other 
portions of the fan assembly. First, the detached blades impacted the engine casing and began to interact 
with other nearby blades. Fragments from impacted blades became stuck between the other blades and the 
fan containment, causing further plastic deformation of the casing. 

While the fan assembly rotated through one revolution following the impact, the detached blades caused 
the fan assembly to move off center, pushing blades on the opposite side of the impact zone into the fan 
casing, resulting in high plastic strains and disintegration. 

The response with no casing was much more uniform in nature. While time and computational resources 
prevented the simulation from running until a quasi-steady state was reached, it was clear from the no 
casing scenario that the response of the system was smooth and approaching a limit cycle as the center of 
mass shifted. 
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Fig. 9. Fan blade tips impacting and sliding against the casing across from the impact zone. 

In contrast, the simulation with the casing present shows the effect that blade -casing contact imposes on 
the vibration response. The fan assembly experienced quick changes in its motion as the complex 
interaction between blades and the casing occurred. Progressive damage of the fan blades and other 
components is shown at various times throughout the analysis in Fig. 10. 
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Fig. 10. Progressive damage at selected times. 


The dynamic oscillation response of the fan assembly with the casing surrounding the fan assembly 
was compared with the simulation without the casing in order to observe the effect that the casing had on 
the vibration. The x- and y-displacements of the center of the hub were plotted through time for each 
model, generating two orbit plots shown in X-Displacement (m) 

Fig. 11. 




X-Displacement (m) 

Fig. 11. Center hub displacement orbit plots. 
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As the unbalanced rotation continued, addition blades around the entire fan assembly made contact with 
the fan casing and experienced yielding and failure at the blade tips as well as along their length. The 
casing was heavily deformed due to this interaction, causing material failure in several locations as a 
result of high strain levels (Fig. 12). 



Fig. 12. Casing plastic deformation after impact. 

Additional simulations with the casing present were conducted in which the number of elements used to 
represent each fan blade was varied. This established the mesh dependence in fan blade damage 
assessment as well as the overall predicted damage. Fig. 13 shows the results of three simulations 
utilizing different fan blade mesh densities immediately after the initial hit. 


(a) 
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The low mesh density model predieted that the initial hit only eaused detaehment of one blade, while the 
medium and high density models predieted detaehment of two blades. This eomparison indieated that a 
degree of fineness equal to or greater than the medium density model was neeessary throughout the 
middle seetion of the blade to aeeurately eharaeterize the bird/blade interaetion between the meshless 
partieles and the struetured Lagrangian mesh used to eonstruet the fan blades. Fig. 14 shows the results of 
the three simulations further in the analysis after several rotations of the fan assembly. 


9 Detached Blades 


15 Detached Blades 




Low density Medium density 

^ ( 1 2420 elements) ( 1 6524 elements) 

18 Detached Blades 



c) High density 
(21600 elements) 

Fig. 14. Damage after several rotations for (a) 12420, (b) 16524, and (c) 21600 blade elements. 


The coarse density blade mesh model predicted a lower rate of damage than the other models, indicating 
detachment of about half of the blades as predicted by the medium and high mesh density models at the 
same time step. The high density model determined a few more detached blades than the medium density 
model. All models predicted substantial damage to the blade tips as they sheared against the casing. 

The larger elements of the coarse mesh model failed at the blade tips as the simulation progressed, 
resulting in a large gap between the remaining blade material and the casing. This gap prevented further 
interaction between the blades and casing, which particularly lead to prediction of less overall damage 
and fewer detached blades. The medium and high mesh density models experienced smaller gaps as 
elements at the blade tips failed, allowing further interaction between the blades and casing. Therefore, a 
level of mesh refinement equal to or greater than that of the medium density model was deemed necessary 
at the blade tips to accurately characterize damage in the entire forward fan section. 

A preliminary assessment of the damage predicted by the model with a real turbofan engine subjected to 
bird strike is shown in Fig. 15. This comparison supported the argument that the widespread damage to 
the fan stage in its entirety, as determined by the numerical model, is possible in cases of real bird strike 
damage. 
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Fig. 15. Comparison of predicted damage with a real bird strike case. 

7 Hybrid Fan Blade Modeling Approach 

In modem engine designs, eomposite fan blades are manufaetured using unidireetional earbon fiber- 
reinforeed epoxy pre-preg material [14]. Layers of the material are arranged and built up on preforms and 
then bonded together to ereate the fan blade stmeture. Metal reinforeements, sueh as a titanium sheath 
that wraps around the leading edge of a eomposite blade, are added to improve the impaet resistanee. The 
resulting assembly is referred to as a hybrid stmeture. 

7.1 Composite Blade Damage 

Failure in hybrid stmetures is eharaeterized by numerous eomplex damage meehanisms oeeurring over 
multiple interaeting length seales and among materials with vastly dissimilar properties. Material models 
that provide a ply-level assessment of eomposite damage were explored in this initial investigation. 
Models that attempt to prediet miero-level fiber and matrix failure, ealled mieromeehanies models, are the 
subjeet of eurrent researeh and will be explored in future studies. 

Ply-level failure meehanisms are broadly defined in terms of strength in the longitudinal (fiber), 
transverse (in-plane perpendieular to fiber), and through-thiekness (out-of-plane perpendieular to fiber) 
direetions, ineluding shear properties relative to these direetions. Failure of an individual ply layer that 
degrades the stmetural performanee but does not eause a eomplete loss of load-earrying eapaeity is 
referred to as damage. Onee all of the ply layers in a given loeation are degraded, the eomplete failure 
oeeurs. The latter is referred to as failure. 

An orthotropie material model with a strength-based failure eriterion was used to model the behavior of 
the eomposite portion of eaeh blade. Ply failure was determined by ealeulating generalized stresses in 
both the longitudinal and transverse direetions for eaeh layer. Onee the stress state satisfied the failure 
eriterion, the load-earrying eapaeity of the layer was removed from the ealeulation, thus eapturing damage 
through the redueed stiffness and strength properties. 

This material was implemented within shell elements with a laminated plate formulation. Integration 
points in the through-thiekness direetion represented eaeh ply layer and allowed eaeh layer’s eontribution 
to the laminate stiffness to be ineluded via the Labotto integration rule [13]. Onee eaeh layer satisfied the 
failure eriterion, the element was removed, thus underlining the loeal eomposite failure. 

Material properties for IM7/8551-7 graphite fiber eomposite were used for the material definition. A 
eomposite layup of [0°/+45°/-45°/0°] was used for this initial study, with the 0° fiber orientation direeted 
along the length of the blade. The thiekness varied aeross the blade from 6.35 mm (0.25 in) at the tip and 
50.8 mm (1 in) at the base, whieh is eonsistent with hybrid fan blade dimensions [14]. 
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7.2 Leading Edge Delamination 


Another critical damage mode for hybrid structures is ductile failure and separation of the metal leading 
edge from the composite edge when subjected to impact loads. In order to capture the damage to the 
leading edge and separation of the sub-structure from the composite blade, the edge was explicitly 
modeled using shell elements that wrapped around the blade edge. 




Fig. 16. Hybrid fan blade model with leading edge explicitly modeled. 


The decohesion technique was chosen for both bonding the leading edge to the blade and capture de- 
bonding of the component. A mixed-mode bilinear traction-separation law was explored that reproduced 
softening of the interlaminar bond between surfaces as they displace, followed by failure as they separate 
a critical distance. 


In this model, the displacement between bonded shell elements at which the softening began was 
determined by the equation 


^8 _ I 


5 '“ = _ 5 '“ = — 


( 11 ) 


where S'** is the displacement for onset of softening, and are the mode I and mode II damage 
initiation separations, ^ is the mode mixity, T and S are the peak normal and tangential tractions, and 
and E^ are the cohesive normal and tangential stiffness, respectively. The displacement at which failure 
occurred was defined by the equation 




1 



( 12 ) 


where is the ultimate displacement, ^IC and are the mode I and II energy release rates, and ct is 
the mixed-mode exponent. 
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This cohesive model was implemented in the explicit FE code using tied a constraint between the titanium 
leading edge and the portion of the composite surrounded by the metal component. With the tied 
constraint approach, parts being coupled do not need to share similar geometries or meshes. This provides 
a convenient method of defining cohesive behavior between parts with very different shapes. 

Cohesive parameters for the composite/metal bond were not available at the time of this study. As an 
approximation of the peak tractions and the mode I and II energy release rates, properties from carbon- 
fiber composite delamination experimental studies were obtained from the literature [15] and decreased 
by one-fourth before being applied, with the assumption that the composite/metal bond would form a 
weaker interface. 

The titanium leading edge incorporated an elastic-plastic material model with a strain-based failure 
criterion. Properties for Ti-A16-4V ASTM Grade 5 titanium alloy were used to define the leading edge 
behavior. 

7.3 Pre-stress Analysis 

Centripetal loads were applied to the hybrid fan blade using the same explicit DR analysis approach 
utilized for metal blades. Again, it was important to decrease the convergence tolerance as much as 
possible to arrive at a pre-stressed state sufficiently close to the true stresses for a given rotational speed 
to reduced unwanted stress oscillations at the beginning of the transient analysis. 

8 Hybrid Fan Impact Modeling 

The fan assembly was formed using an array of hybrid fan blades arranged around a central hub and 
connected via rigid constraints at their bases. Half of the fan assembly was modeled with a total of 9 
blades. The 4 lb. bird model developed earlier was introduced to impact at about the mid-span of the 
blades. Contact was defined between the particles and the hybrid blades using the same penalty-based 
approach. 

The fan assembly was set rotating with an initial speed of 2400 rpm, a typical maximum speed for similar 
high bypass ratio turbofans with wide-chord fan blades. The bird was assigned an initial relative velocity 
of 100 m/s. 



Fig. 17. Explicit FE model of a bird strike scenario for a hybrid fan assembly. 
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8.1 Bird Impact Results 


The bird was sliced into two parts upon impact with the first blade. Part of the bird material remained in 
contact with the blade, while the remaining material continued to impact additional structures. The soft 
body was separated into multiple impactors due to slicing by subsequent blades, with portions coming 
into direct contact with four fan blades. Contact between the structures and the bird involved both a 
slicing action with the blade edge and a blunt impact with the rear face of the blade. The impact event is 
shown in Fig. 18. 

None of the hybrid blades experienced failure in any part of the structure as a result of the impact. 
However, critical damage mechanisms were captured throughout several of the blades impacted by the 
bird. The damage analysis showed fiber tensile failure occurring at the base of two blades as a result of 
high stress at their respective roots, as shown in Fig. 19. 



Fig. 18. von Mises stresses induced by the bird impact at progressive time steps. 


These blades were subject to impact with the largest portions of the bird. The plot shows the percentage 
of plies intact/failed, with blue indicating 100% of plies intact and red indicating 100% of plies failed. 
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Fig. 19. Fiber tensile failure initiated at the root of two impacted blades. 

In addition, matrix failure was predicted throughout each of these blades, especially where the bird 
impacted the face of the blade rather than the leading edge. This damage is shown in Fig. 20. 



It was observed that damage was more prevalent in the outer plies. The outer plies are subjected to higher 
normal strains and stresses when subjected to bending caused by the bird impact. This tendency towards 
outer ply damage is shown in Fig. 21. 
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Fig. 21. Matrix and fiber tensile damage within each layer for an impacted hybrid fan blade. 


Leading edge delamination was predieted in several blades by the eohesive model. The delaminated 
region is shown as red in Fig. 22. Blue indieates fully bonded areas. Other regions indieate softening prior 
to full de -bonding. 
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Fig. 22. Leading edge delamination following bird impact. 


9 Comprehensive Crashworthiness Methodology 


This work and complementary studies have demonstrated that interactions among forward section 
components contribute to additional damage, not sufficiently captured by a stand-alone assessment of bird 
impact on a structure [16]. As an additional step towards developing a comprehensive crashworthiness 


NASA/TM— 2014-2 1 8397 


21 



methodology for hybrid engines, a fully eoupled fluid-struetural module is proposed. The requirements 
for this module and the platforms utilized will be detailed in the next report. 

10 Conclusion 

Researeh efforts foeusing on the development of a detailed predietive modeling methodology for 
eomprehensive analysis of bird strike and post-impaet damage assessment in struetural -hybrid jet engine 
forward seetions are reported. The development of a transformative, aerospaee propulsion design platform 
for evaluating the erashworthiness of future advaneed struetural engines is one of the main objeetives of 
these studies. The novel methodology proposed offers the ability to eapture the extent of intrieate soft 
foreign objeet damage to fan blades resulting from direet initial impaet, subsequent eollisions with a 
number of other engine forward seetion parts and sub-eomponents as well as assoeiated fragmentation 
eontainment analysis. The full dynamie response prognosis in this methodology further eneompasses 
blade-out, fan assembly unbalaneing and vibrations, and subsequent engine easing shearing and shaving 
damage analyses, whieh remain unaddressed in many previous studies. 

A novel numerieal method for predieting damage and failure in hybrid jet engine fan blades subjeet to 
soft impaet has also been developed. Critieal damage meehanisms have been eaptured, sueh as eomposite 
fiber and matrix failure meehanisms, leading edge delamination, and material failure. This work shows 
promise for extending the forward seetion engine damage methodology to inelude hybrid fan struetures. 

The outlook of the investigations will inelude validations of the modules developed thus far by means of 
eomparison with data obtained from variety of erashworthiness experiments. Overall, it is aimed to 
develop a eomprehensive analysis platform for the assessment of hybrid engine designs that will reduee 
relianee on physieal tests and ean support future developments of advaneed aerospaee propulsion systems. 
Fluid-struetural eoupling will further enable the methodology to assess the performanee loss eaused by 
struetural deformation and failure. 
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